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Using a density functional approach based on a Skyrme interaction, thermodynamic properties of 
finite nuclei are investigated at non-zero temperture. The role of a momentum dependent isovector 
term is now studied besides volume, symmetry, surface and Coulomb effects. Various features 
associated with both mechanical and chemical instability and the liquid-gas coexistence curve are 
sensitive to the Skyrme interaction. The separated effects of the isoscalar term and the isovector 
term of momentum dependent interaction are studied for a modified SKM(m* — m) interaction. 
On ' The frequently used Skyrme interaction SLy4 is one of the cases considered and is shown to have 

' better features for neutron star studies due to a larger symmetry energy. 
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■ I. INTRODUCTION 

' The study of nuclear properties at finite temperature is important for nuclear astrophysics, medium energy heavy 
ion collisions and future experiments at the Rare Isotope Beam Facility (RIB). It is also of general interest because 
^SJ ' of its relation to the study of strongly correlated fermions and to phase transitions in multi-component systems, and 
in particular to these studies in two component systems. Density functional theory has been extensively used as an 
approach to such studies in many areas of many body physics. In nuclear physics density functional theory based on 
a Skyrme interaction approach has been frequently used. Here such an approach will be employed in our study of 
J-{ nuclear properties at finite non-zero temperature. Most applications of density functional theory are for systems at 
O T = 0. 

The nuclear system has terms that set it apart from other two component binary systems. These include Coulomb 
terms from the charged proton component, symmetry terms from the nuclear symmetry energy, surface energy terms 
and momentum dependent terms from a velocity dependence of the nuclear interaction. Questions related to nuclear 

■ stability should involve these terms since it is their interplay that determines important features. For example, fission 
processes are in part an interplay of surface versus Coulomb interactions. The valley of nuclear stability in N and 
Z involves an interplay of symmetry versus Coulomb interactions. While the Coulomb energy favors systems with 
large neutron excess in heavy nuclei, symmetry terms favor systems with equal number of protons and neutrons. 
Nuclei far from the valley of nuclear stability up to neutron drip and proton drip lines are of interest in future 
RIB experiments. Moreover the Coulomb force is long range and the nuclear force is short range and nearly charge 
independent. A charge independent nuclear interaction has an associated isospin sysmmetry. The presence of such 

' terms in the interaction of nucleons make the nuclear case a unique binary system to study. Moreover, realistic 
, nuclei are finite and contain surface terms so that phase transitions in finite systems can be studied. In this paper 
IIV ' we incorporate an additional isovector component in the momentum dependence of the interaction, extending our 
. ^ [ previous studies [ESS] of properties of heated nuclei. In our previous work we considered an isoscalar momentum 
dependent term besides volume, surface, non momentum dependent symmetry terms and Coulomb interactions. A 
momentum dependent isovector term and related isovector effective mass effects are important in a study of properties 
associated with symmetry terms between protons and neutrons. How various terms in nuclear two component systems 
manifest themselves at non-zero temperature and how they affect the equation of state, compressibility, chemical and 
mechanical instability and associated liquid gas coexistence curve are investigated. Without Coulomb terms, the 
coexistence curve and chemical and mechanical instability curves are symmetric about proton fraction y = 1/2. 
However a stable nucleus such as Pb^°^ at zero temperature has proton fraction y = 82/208 « 80/200 = 0.4 far from 
being symmetric in proton number 82 and neutron number 126. When Coulomb terms are included a large asymmetry 
appears in curves associated with finite temperature properties. The symmetry terms partially restores the symmetry 
of the curves. We also extend our previous calculation by considering a commonly used Skyrme interaction labeled 
SLy4 and discuss its properties. 

The phase diagram of a one component system is a simple curve of pressure versus density determined by a 
Maxwell construction. The phase structure in a binary system is considerably more complex than one component 
systems because of an extra degree of freedom associated with the additional component. The binodal surface is now 
determined by pressure, temperature and proton fraction. The proton fraction can be different in the denser liquid 
phase than the less dense gas phase because of the short range of the nuclear symmetry energy. In a nuclear two 
component system, isospin fractionation [J, [^, @, 0j B is an example where the monomer gas has a large neutron 
excess. Ref.Q is a recent reference that contains further references to isospin fractionation. The liquid gas phase 
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transition was first treated as a one component system jio'l and later extended to two components [Tl| using a Skyrme 
interaction. Ref.[T^ used a relativistic mean field model and studied the role of the symmetry energy on the phase 
structure in detail. Inclusion of Coulomb and surface terms can be found in Refs. [H, Q, Q where large asymmetries 
were shown to appear from the Coulomb term. Pawlowski 13] has also considered the role of the Coulomb term. 
Some further studies of the nuclear phase transition can be found in Refs.[i3,[il[ii,[i3,[ii,[ii,[2S[2l|,[2l[2l[23,i^ 
Ref. [i^l also includes a discussion of momentum dependent terms, both isoscalar and isovector. Our results will differ 
somewhat from those of Ref. ^2§\ because we consider a Coulomb term which leads to asymmetries in various quantities 
in proton fraction. Without a Coulomb interaction symmetry exists around a proton fraction y — 1/2 for nuclear 
stability and for the coexistence loops. The importance of nuclear isospin symmetry in various features of heavy ion 
collisions can be found in Refs.0, Hg]. 

In Sect. II, we give the main equations that are necessary for understanding the thermodynamic properties of two 
component systems with momentum dependent interaction at non-zero temperature. Sect. Ill is an application of 
the results of Sect. II to questions associated with the mechanical and chemical instability of hot nuclear matter and 
the associated liquid-gas phase coexistence curve. A summary and conclusions are given in Sect. IV. 



II. THEMODYNAIC PROPERTIES OF NUCLEI IN A SKYRME DENSITY FUNCTIONAL 

DESCRIPTION 



A. General Results 



In this section we present results for the thermodynamic properties of nuclear matter which are extended from 
the results of Ref.[l|, 0, Q to now include a isovector velocity or momentum dependent interaction. For a Skyrme 
interaction, this momentum dependent isovector term has some important consequences for both nuclear stability and 
for phase transition properties. 

For a nuclear system of proton {pp) and neutron this gives the local potential energy density as 
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yi/'P P^^^ with 77 = 2/3 when we approximate the Coulomb 
and surface effects as coming from a finite uniform sphere of radius R = r^A^^^ with total charge Z which has the 
Coulomb energy density of Uc = f ^^f^ as discussed in Ref.[3|. We have used the value of i? = 6 fm and AnrQa — 20.0 
MeV. The values for the force parameters used here are given in Table [H We define an effective mass m* as 
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where p ^ pp + pn and y — Pp/p. The upper and lower signs are for proton and neutron repectively. When 
Xl = X2 ^ ^1/2, the m/m* is a pure isoscalar since the coefficient in the term of pq vanishes. The second form 
of Eq.(I21) shows that the effective masses of neutron and proton are the same for a symmetric nuclear system. We 
can also see that the neutron effective mass becomes heavier and the proton mass becomes lighter in a neutron rich 
system while proton effective mass becomes heavier and neutron mass lighter in a proton rich system. 

In Table m the column involving no isovector term is obtained from the column SKM(m* = m) by simply setting 
Xl = —1/2 and X2 = —1/2 and the column with a momentum independent interaction is obtained by setting 
ti — t2 ~ 0. Both these columns leave the other parameters the same. From the table we see that simply changing 
these parameters changes various quantities given at the bottom of the table starting with the effective mass and 
ending with the compressibility. It should be noted that nuclear matter properties for the momentum independent 
interaction and the SKM(to* = m) are the same except for the symmetry energy. Specifically, when two interactions 
have the same effective mass then the energy and pressure are the same in symmetric systems with the potential 
energy of the form of Eq.([T]) (see Eqs.(42), (45), and (46) of Ref. Q). And thus the saturation properties will be the 



3 



TABLE I; Skyrme parameters are in MeV and fm units. First three columns are from the SKM(m* — m) parameter set [H. [27l| 
except xq — —1/6 and X3 — —1/2. The last column is from the SLy4 parameter set [2^. The nuclear matter properties at the 
saturation are also shown. The effective mass is for a saturated symmetric nuclear matter. 





SKM(m* = m) 


No Isovector 


Momentum indep. 


SLv4 




-1089.0 


-1089.0 


-1089.0 


-2488.91 


Xo 


-1/6 


-1/6 


-1/6 


0.834 


I'd 


17270 


17270 


17270 


13777.0 




-1/2 


-1/2 


-1/2 


1.354 


Q! 


1 


1 


1 


1/6 


t'i 


251.11 


251.11 





486.82 







-1/2 





-0.344 


t1 


-150.66 


-150.66 





-546.39 


X2 





-1/2 





-1.000 


Effective mass m* /m 


0.999987 


0.894430 


1 


0.694658 


Binding energy Eb/A 


15.8173 


13.3250 


15.8176 


15.9722 


Fermi energy Ef 


34.5186 


32.0728 


34.5188 


36.7743 


Saturation density po 


0.14509 


0.12994 


0.14509 


0.15954 


Symmetry energy Sv 


18.6080 


23.7451 


24.6730 


32.0038 


Compresibility k 


367.556 


312.281 


367.562 


229.901 



same. The only difference between these two columns is then the symmetry energy. By contrast, by simply changing 
the effective mass, the energy and pressure in the system are changed and thus also the saturation properties. This is 
seen in Table|T]by comparing the SKM(m* = m) column with the column of no isovector interaction. Our purpose of 
just varying xi and X2 or ti and ^2 without varying other parameters was to isolate the specific role of the isovector 
and isoscalar terms in SKM(m* = m) and thus study their role. We see that various properties of nuclear matter 
listed in the table are changed and the degree to which they are changed is presented in the results. Changing xi 
and X2 (or ti and ^2) and also changing the other parameters to fit the saturation properties gives rise to yet another 
Skyrme interaction. Comparisons with many Skyrme interactions can also be made which we hope to do in the future. 

For a nuclear system with protons and neutrons with the interaction given by Eq. ([1]) , the non-degenerate Fermi gas 
limit [U, 0, Q leads to the following set of equations. The chemical potential has a behavior determined by 
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The equation of state has a form given by 
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The effective mass m* and thus are, in general, isospin dependent |26l . |2E 



For 77i/m* = l + ogp = 1 + ap + bpg = 1 + p(a + 6[l± (2y- l)]/2) with pq = p[(l=Fl)/2±y] and A, = J2Trhym*T 
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The ap term in the density dependent effective mass comes from the isoscalar term of momentum dependent interaction 
and the bpq term comes from the isovector term of momentum dependent interaction. Using a and 6, we can study 
the variation of P and pq with p and y. We can study the behavior of thermodynamic quantities at a fixed P using 
dP = from Eq.®, 
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This equation gives the condition determining ym{p) where both dP/dy = and dp/dy — 0, 
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The Urnip) is the value of y with lowest pressure P for a given density p and temperature T and is indepenent of p 
for a momentum independent Skyrme interaction with = —1/2 and /? = as considered in Ref.[li, i^]. The term 
and the density dependent effective mass for a momentum dependent Skyrm force introduce a small p-dependence in 
Urn- The isovector momentum dependent term {xi and X2 terms) introduce large p-dependence in ?/„ as we will see 
later. The ym{p) curve cross the coexistence loop at the equal concentration point ys where the liguid and gas phases 
have the same y value. 
From Eq.®, 
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This equation determines the various curves of dpq/dp ~ 0, dpq/dy = 0, chemical instability curve, etc. 
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FIG. 1: Equation of state P(p) for fixed y = Q, 0.2, ym{p) (the proton concentration with minimum pressure defined by 
Ea. (llO|l l for thick lines from top to bottom and y = 0.5, 0.8, 1.0 for thin lines from bottom to top. The solid line is for the 
modified SKM(m* = m) interaction with both isoscalar and isovector momentum dependent terms, the dash-dotted line is for 
a momentum dependent but isovector independent term, and the dashed line is for momentum independent terms only. The 
lowest thin dashed curve overlaps here with the lowest thin solid curve which have y = 0.5. 



III. APPLICATIONS TO NUCLEAR MECHANICAL AND CHEMICAL INSTABILITY AND THE 

LIQUID-GAS PHASE TRANSITION 

A. Equation of State with Isoscalar and Isovector Momentum Dependent Terms 

The equation of states {P{p) curve) for various proton fraction y at a fixed temperature of T — 10 MeV are shown 
in Fig[T] The solid lines are P{p) for momentum dependent interaction (modified SKM(m* = m)) with both isoscalar 
and isovector terms. The dash-dotted lines are for a momentum dependent Skyrme interaction with isoscalar term 
only and the dashed lines are for a Skyrme interaction without momentum dependent terms which are the same 
cases considered in Ref.Q. For various values of pressure P at T = 10 MeV, the curve of y{p), fJ-q{y), and Pqip) 
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are shown in Fig[2l FiglSl and Fig[4] respectively. The basic behaviors of the curves are the same independent of 
which force parameters are used. For example, the y{p) and fJ-q{p) curves of Figsl2]and|3]for low pressure P exhibit 
a closed loop separated from a vertical line at low density for all three cases. Since the overall behaviors which are 
indepenent of force parameters are studied in previous papers [H, 0) Q i the discussions below are concentrated more 
on the differences between different force parameters. 

The pressure for a given density p and y increases from the corresponding value for the case of a momentum 
independent Skyrme interaction (dashed curve) by the momentum dependent isoscalar term (dash-dotted curve) and 
decreases much more by the isovector term (solid curve). In Fig.l, the crossing of curves for different y values appears 
explicitely for y = 0.2 and 0.5 with a momentum dependent isovector term (solid lines). This crossing is related to 
the p dependence of y™, the proton concentration with minimum pressure defined by Eg. ljlOp . which comes from the 
momentum dependent terms and the three body isovector term. 

For y = and with an isovector momentum dependent term the P{p) curve has a spinodal instability region in the 
density range of about 1/5 < p/po < 1/3 of nuclear saturation density po and around P = 0.08 (0.0755 ~ 0.0834) 
MeV/fm^. In Figs E] and H the curves for P = 0.08 MeV/fm^ show a gap in the range of 0.0401 < p < 0.0598 fm^^ 
due to the spinodal behavior for y = 0. From solid curves of Fig[l]we can see the horizontal cut of P = 0.08 McV/fm"^ 
gives negative value of y in this density range. This spinodal instability of P{p) curve for y = causes the cut of 
coexistence curve by the y = line which can be seen in Figs [5] and [6] 

For y = 1/2 the P{p) curve for momentum dependent interaction with both isoscalar and isovector terms (solid 
curve) overlaps with the curve for momentum independent interaction (dahed curve). The SKM(m* = m) interaction 
used here (Table U), with both isoscalar and isovector momentum dependent terms, has the effective mass of m* = ni 
for y = 1/2 which is the same value with the momentum independent interaction case. The y{p) curves at fixed P 
in Fig[2] show that the curves for a momentum independent force (dashed line) are tangent to the curves for fully 
momentum dependent force (solid line) at y = 1/2. Similary the Pq{y) curves in Fig[3]show that the dashed curve 
crosses the corresponding solid curve a,t y — 1/2. In Fig|4] this crossing appears at the point with large p which 
corresponds to y = 1/2. Since m* = m at y = 1/2 for a full momentum dependent force the P{p) curve aX y = ym{p) 
is closer to the curve for momentum independent force than to the curve for a momentum dependent interaction 
without an isovector term. 

Figsl2] and [6] each show that the momentum dependent isovector term introduces a large p dependence in the ym 
of Eq. (|10p . the proton concentration y having minimum pressure for a given density p. The ym{p) is independent 
of the density p for a momentum independent Skryme force. The isoscalar momentum dependent term introduces a 
small increase of ym with increasing p while the isovector momentum dependent term introduces a large decrease of 
ym with increasing p. 

FigsH] - ini show that the basic behavior of various curves remains qualitatively the same for the various cases 
considered here. However the isoscalar momentum dependent term and the isovector momentum dependent term 
move in an opposite direction for the various curves compared to the curves of a momentum independent case. The 
isovector term has a larger effect than an isoscalar term. 

The y{p) curve shown in Fig[2]is the horizontal cut of a fixed P in FiglT] Fig[2]shows the closed loops of y{p) for low 
pressure P is reduced to a much smaller loop by the isoscalar momentum dependent term (dash-dotted line) from the 
momentum independent interaction (dashed line). The isovector momentum dependent term makes the loop much 
larger. Thus, the full momentum dependent interaction (solid line) enlarges the loop from the momentum independent 
loop. For higher P which has no closed loop, the isoscalar momentum dependent term makes the density p for a 
given y and P smaller than p for momentum independent force and the isovector momentum dependent term makes p 
larger. The full momentum dependent term makes p for a given y and P larger than the p for momentum independent 
force. The proton fraction y for a given p and P is farther away from ym for full momentum dependent interaction 
than other interactions due to the smaller value of symmetry energy. The solid curve for P — 0.08 MeV/fm'^ shows a 
cut off of the curve by y = line for the density range of 0.04 0.06 fm~'^. 

Fig[3] and |4] show that the chemical potential has a somewhat complicated relative effect arising from momentum 
dependent terms. Except in the regions with medium y in Fig[3]or larger value of p in Fig[?l the isoscalar momentum 
dependent term makes the chemical potential higher while isovector momentum dependent term makes it lower. The 
closed loop for low P in Fig|3] is reduced to a much smaller size by momentum dependent isoscalar term (dash- 
dotted curve) from a loop for a momentum independent interaction (dashed curve) while it is enlarged by momentum 
dependent isovector term (solid curve). The low pq portion of the curves for higher P in Figdl which corresponds to 
the region with lower y value, shows the density for a given pq becomming much smaller by the isoscalar momentum 
dependent term from the value for the momentum independent case while it becomes much larger by the isovector 
momentum dependent term. The horizontal straight thick solid line (the full momentum dependent SKM(m* = m) 
force) at low pp for P = 0.08 MeV/fm"^ in Fig 0] represents the cut due to the condition of y > for the density range 
of 0.04 < p< 0.06 fm-3. 
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Density, p [fm ^] 

FIG. 2: The y{p) dependence for fixed pressure of P = 0, 0.015, 0.05, 0.08, 0.2, 0.3, 0.4, 0.5 MeV-fm"^. The curves are the 
same as in Fig.l. The thin horizontal lines are for ym{p) curve define by Ea. (|10|l . 

B. Liquid-Gas Phase Transition and the Coexistence Curve and Instability 

For a one component system the coexistence curve is a line obtained by the famihar Maxwell construction as 
already noted. For a two component system the coexistence region is a surface obtained as follows. The condition 
for coexistence between the two phases requires the proton chemical potentials to be the same in two phases and, 
similarly, the neutron chemical potentials must be the same in the two phases at a given pressure and temperature. 
Note that the proton fraction or neutron fraction need not be the same in each of the two phases. The two phases 
have the same proton fraction at the line of equal concentration ije which is the line of intersection between the 
coexistence surface and the surface of j/m satisfying Ea. (fTO| . In fact, the liquid phase should be a more symmetric 
system than the gas phase because of the symmetry potential as seen in Refs.[ll, This observation goes under the 
name of isospin fractionation. 

Figs l5] and [5] show features of the coexistence curves. Also shown in FiglHlare the mechanical and chemical instability 
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Proton Fraction, y 

FIG. 3: Chemical potential fip (thick lines) and fj,n (thin lines) as a function of y for P = 0.015, 0.05, 0.08, 0.2, and 0.5. The 
curves are the same as in Fig.l 



loops and curves of dfiq/dp — 0. In Figs (5] and [5] the small kink at lowest y side of the coexistence loop for the 
interaction with an isovector momentum dependent term might have come from a numerical problem. For SLy4 
(curve with + sign in FigO, there is a small island of coexistence near the yE region. It is not clear if this is a 
numerical problem or a physical feature. Further detailed studies are required. 

FigUl shows that the coexistence loops in the neutron rich side (smaller y side) are smaller than the loops in the 
proton rich side for three cases of SKM(m* — m) parameter sets. For the SLy4 parameter set, the coexistence loop 
at the neutron rich side is larger than the loop at the proton rich side. Ref. Q shows that the Coulomb interaction 
moves yE to a. smaller proton fraction y and makes the coexistence loops smaller in the neutron rich side than in the 
proton rich side. The SKM(m* = m) parameter sets used in Ref. 2] has no isovector three body term {x-^ = —1/2) 
while SLy4 parameter sets has a negative three body isovector term. 

In Fig[5]we can see the instability curves are inside of the coexitence loop for all three force parameter sets. The 
mechanical instability loops are inside of the chemical instaility loops and tangent to each other at the point of 
lines (the horizontal line in this figure). The chemical instability loops are inside of the coexistence loops and tangent 
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Density, p [fm "^J 

FIG. 4: Same as Fig. 3 but as a function of p. The straight horizontal portion of the thick sohd curve (momentum dependent 
force with both isoscalar and isovector terms) for P = 0.08 MeV/fm"^ at low /ip just represents the cut in the density range of 
0.04 < p < 0.06 fm"^ due to the condition of y > 0. 

to each other at the critical point, i.e., the point with highest pressure on the coexistence curve. 

In Figs [5] and [6l the loops for the isovector momentum dependent case (solid line) has a cut at y = which is 
related with the spinodal behavior of P{p) curve for y — Q in FiglT] The smaller symmetry energy allows the system 
to be more asymmetric and thus to reach the boundary of y = or y = 1. The SKM(m* = m) force with both 
isoscalar and isovector momentum dependent terms has the smallest symmetry energy while the SLy4 force has the 
largest symmetry energy among forces used here (see Table H]). FigO shows the coexistence loops for SLy4 are more 
symmetric than the loops for SKM(m* = m) parameter sets and the full SKM(m* = m) force has the most asymmetric 
coexistence loops. The cut of coexistence loop by j/ = indicates that the SKM(m* = m) force with both isoscalar 
and isovector terms has too small of a symmetry energy for a description of nuclear system even though we have used 

= —1/6 instead of a;o = 0.5 to have a larger symmetry energy. 

The pressure of the coexistence curve on the neutron rich side for a given value of y (Fig[5]) becomes higher due to 
the momentum dependent isoscalar term (dash-dotted line) than the pressure for momentum independent interaction 
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Proton Fraction, y 

FIG. 5: P{y) plot for coexistence curve. The curves are the same as in Fig.l. The curve with + sign is for SLy4 parameter. 

(dashed line) and also moves the peak to higher y. By contrast the momentum dependent isovector term (solid line) 
makes the pressure lower and moves the peak to lower y. The effects of isoscalar and isovector term work in opposite 
direction for the coexistence curve on the proton rich side compare to the effects for neutron rich side. For SLy4 the 
system without an isovector or isoscalar momentum dependent term does not saturate at the usual normal density 
region. However, the isovetor momentum dependent term in SLy4 has the same sign as the isoscalar momentum 
dependent term thus both isoscalar and isovector terms would make the pressure of the coexistence curve of neutron 
rich side higher than momentum independent case. This might also make the coexistence loop in neutron rich side 
larger than the loop in proton rich side. In Figl6]the loops for the momentum dependent isoscalar case (dash-dotted 
line) is smallest and thus mostly inside among the loops while the loops for the momentum dependent isovector case 
(solid line) are largest and thus mostly outside. The loops for the momentum independent case (dashed line) are 
tangent to the loops for the momentum dependent case with both isoscalar and isovector terms (solid line) at y = 1/2 
where m* — m for both cases. The curves of dfiq/dp — (thin curves in Figl6]) also show that the isovector momentum 
dependent term has an opposite and larger effect than the isoscalar momentum dependent term. 
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Density, p [fm ^] 

FIG. 6: The closed loops are for coexistence curves, chemical instability boundary curves, and mechanical instability curves 
from outmost loops to inside loops at T = 10 MeV. Also shown by thin lines are dfiq/dp — for proton (opened downward) 
and for neutron (opened upward). The horizontal lines are the ym{p) determine by Ea. (|10|) . 



IV. SUMMARY AND CONCLUSIONS 



A density functional theory based on a Skyrme nuclear interaction is used to study the temperature dependent 
properties of a two component system of strongly interacting protons and neutrons. The paper extends our previous 
studies by including an isovector momentum dependent interaction. As before, Coulomb and finite size surface effects 
are also contained in the description of various properties. The nuclear interaction has a velocity or momentum 
dependent term and in a medium this momentum dependence can be discussed as an effective mass to lowest order 
in it. The momentum dependence has both isoscalar and isovector terms. Our prior study in Ref.0] suppressed the 
isovector part by using a particular choice of Skyrme parameters. This greatly simplified the original calculations of the 
stability and coexistence properties. Here we now explore the mechanical and chemical instability of a two component 
system and also the coexistence features associated with a liquid gas phase transition incorporating an isovector 
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momentum dependence besides other terms mentioned above. While quahtative results have a similar behavior with 
and without a momentum dependence (isoscalar and isovector), significant quantitative differences exist. Moreover, 
we extended our calculations to include the frequently used Skyrme interaction called SLy4 [28l |. 

Each term, Coulomb, volume, surface, symmetry, effective mass or momentum dependence, both isoscalar and 
isovector, play a unique role in determining the finite temperature stability properties and coexistence features as 
shown in various figures in this paper. For example, the Coulomb interaction leads to a pronounced asymmetry in 
the coexistence loops shown in FigjSl Moreover the equal concentration point tje, where the liquid and gas have the 
same proton fraction, moves away from tje = 1/2 without Coulomb forces to a point where the proton fraction is close 
to the valley of stability when Coulomb forces are included. For example ue = 0.410 for a momentum independent 
force, 0.414 for a isoscalar momentum dependent force, 0.379 for a isovector momentum dependent force, and 0.402 
for SLy4. The yE with Coulomb included is somewhat insensitive to the nuclear force. However, the pressure at ue is 
more sensitive. Specifically, P — 0.0138 MeV/fm'^ for momentum independent case, 0.0236 for isoscalar case, 0.0127 
for isovector case, and 0.0262 for SLy4. 

A model with no Coulomb terms would have identical loops in y around the point of intersection yE = 1/2 of equal 
concentration. The SLy4 interaction has the interesting property of making the left coexistence loop of lower proton 
fraction more pronounced than the right loop of high proton fraction. This particular feature in the height of the 
right versus left loop is not observed with the other SKM(m* = m) interaction discussed in this paper. In particular, 
the right loop is still more pronounced or higher than the left loop for SKM(m* — to). For SKM(to* = to), the 
isoscalar momentum dependent term raises the pressure of the coexistence loop on the neutron rich side while the 
isovector momentum dependent term, which has an opposite sign from the isoscalar term, lowers the pressure. For 
SLy4, both the isoscalar and isovector terms have the same sign and thus would raise the pressure of the coexistence 
loop on the neutron rich side. The SLy4 has also a non zero effect from an isovector three body term which was zero 
for SKM(to* = to). Besides this feature the SLy4 loop at higher y side is narrower in y for a given T and P than 
the SKM(to* = m) loops. From this observation we can conclude that the coexistence and stability behaviors are 
sensitive to the choice of Skyrme interaction parameters. The SLy4 has a higher symmetry energy than the other 
interactions used in this paper as shown in Table |T1 Moreover a higher left loop means that the liquid and gas phases 
can coexist at higher pressures at a given temperature with a more neutron rich gas and a lesser rich neutron liquid, 
i.e., yo < Vl < He < l/'2 for the left loop. A narrowing of the coexistence loop brings the yc and yL closer together 
thereby reducing the proton fraction difference in the gas and liquid phases. 

The momentum dependent isovector term in SKM(to* = to) force has an opposite and larger effect than the 
momentum dependent isoscalar term. A momentum dependent isovector term makes the symmetry energy smaller 
and thus the coexistence loops and the mechanical and chemical instability loops become more asymmetric. The small 
value of symmetry energy for full momentum dependent interaction makes the y value for a given p and P farther 
away from ym,- The pressure for a given p and y at fixed T becomes larger by the isoscalar momentum dependent 
term while it becomes smaller by the isovector momentum dependent term. 

Symmetry energy terms arising from momentum independent and momentum dependent terms tend to restore 
the symmetry of the loops. The coexistence loops and instability loops for SKM(to* = m) with both isoscalar and 
isovector momentum dependent terms have the smallest symmetry energy. Therefore they are the furthest away from 
yE for this force and thus have a cut due to the y > condition. This feature presents a problem for using an 
interaction with too small of a symmetry term in neutron star studies. The coexistence loops for SLy4, which has 
the largest symmetry energy, are closer to its yE than for the SKM(to* — m) parameter sets. This suggests that the 
large symmetry energy terms of the SLy4 are very important in restoring the isospin symmetry in the system and 
also better in neutron star studies. 
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